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In neutral graphene, two prominent cusps known as Kohn anomalies are found in the phonon 
dispersion of the highest optical phonon at q = T (LO branch) and q = K (TO branch), reflecting a 
significant electron-phonon coupling (EPC) to undoped Dirac electrons. In this work, high-resolution 
electron energy loss spectroscopy is used to measure the phonon dispersion around the T point in 
quasi-freestanding graphene epitaxially grown on Pt(lll). The Kohn anomaly for the LO phonon is 
observed at finite momentum q ~ 2/cf from T, with a shape in excellent agreement with the theory 
and consistent with known values of the EPC and the Fermi level. More strikingly, we also observe 
a Kohn anomaly at the same momentum for the out-of-plane optical phonon (ZO) branch. This 
observation is the first direct evidence of the coupling of the ZO mode with Dirac electrons, which is 
forbidden for freestanding graphene but becomes allowed in the presence of a substrate. Moreover, 
we estimate the EPC to be even greater than that of the LO mode, making graphene on Pt(lll) 
an optimal system to explore the effects of this new coupling in the electronic properties. 


Introduction - Kohn anomalies are kinks in the phonon 
dispersion of a material produced by abrupt variations in 
the screening of atomic vibrations by gapless electrons P- 
[3j. The properties of these anomalies are completely de¬ 
termined by the electron-phonon coupling (EPC) and the 
shape of the Fermi surface, and their observation thus of¬ 
fers a window to study the interplay between electronic 
properties and phonon dynamics. 

As a paradigmatic example, in graphite |UZ| and 
graphene SHE], Kohn anomalies are realized as linear 
cusps in the dispersion of the highest optical phonon 
branches at T (the E 2 g phonon) and at I\ (the A\ 
phonon). The form and position of this cusp is deter¬ 
mined by the Dirac fermion dispersion of the electronic 7r 
bands around K and the Fermi level p, ss 0. These two in¬ 
plane phonon branches are commonly accepted to be the 
only ones with a significant EPC in graphene m- Within 
a tight-binding picture, these are the only branches that 
modify the nearest neighbor hopping integrals, so other 
in-plane phonons require a different mechanism to couple 
to electrons. 

For the case of out-of-plane phonons, the EPC is 
strongly constrained by the presence of mirror symme¬ 
try with respect to the horizontal plane which forbids a 
first order coupling to electrons, at least in free-standing 
graphene. When graphene is supported by a substrate 
a first order EPC becomes allowed in principle p~7l Il8j . 
but it has remained unknown whether its strength is large 
enough to produce any significant effect. In particular, 
for the out-of-plane optical mode (the ZO mode), it has 
been speculated that a finite coupling could be respon¬ 
sible for a Peierls instability to a spontaneous buckling 
of the lattice mm, which has not been observed so 
far. A coupling of the ZO mode to electrons would in¬ 
fluence many physical properties, since the ZO has sig¬ 


nificantly lower energy than the in-plane branches with 
strong EPC. The observation of a Kohn anomaly for 
the ZO phonon [2T] would represent a definitive proof 
of the existence of the EPC for this phonon branch and 
would allow an estimate of it. Herein, we demonstrate 
the existence of precisely this Kohn anomaly by means 
of high-resolution electron energy loss spectroscopy [22] 
(HREELS) measurements on monolayer graphene grown 
on Pt(lll). In addition, we also present evidence of the 
expected Kohn anomaly for the E^g phonon branch. 

Graphene on Pt(lll) is an ideal playground for inves¬ 
tigating the possible emergence of a ZO Kohn anomaly, 
as it is characterized by the weakest graphene-substrate 
interaction [23]. The graphene-Pt distance (3.31 A) lies 
close to the c-axis spacing in graphite. Moreover, the 
electronic structure of graphene on Pt(lll) resembles 
that of isolated graphene [24] > with the Dirac fermion 
dispersion of 7r bands preserved. Angle-resolved pho- 
toemission spectroscopy (ARPES) experiments [24] do 
not show any significant hybridization of the graphene 7r 
states with Pt d-states, which simply superpose in energy 
with minimal interaction between them. This is in con¬ 
trast to the case of, for example, graphene on Ni(lll), 
where the hybridization of the graphene 7r-states with the 
Ni d bands has a very strong effect on the ir bands [25] . 
The graphene sheet grown on Pt(lll) is p-doped with a 
shift by 0.30 ± 0.15 eV of the Fermi level from the Dirac 
point. 

Experimental methods - Experiments were carried out 
in an ultra-high vacuum chamber operating at a base 
pressure of 5 • 10~ 9 Pa. The sample was a single crystal 
of Pt(lll). The substrate was cleaned by repeated cy¬ 
cles of ion sputtering and annealing at 1300 K. Surface 
cleanliness and order were checked using Auger electron 
spectroscopy and low-energy electron diffraction mea- 
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FIG. 1. HREELS intensity plot for phonon dispersion of 
graphene/Pt(lll). Phonon modes have been recorded in 
HREELS spectra acquired in off- specular geometry with the 
sample oriented along the P — M direction. The incidence 
angle is 0; = 80° with respect to the surface normal and the 
primary electron beam energy is E p = 20 eV. To put in ev¬ 
idence only inelastic losses due to phonons, the tail of the 
elastic peak is grayed out for clarity. The LO, LA, ZO and 
ZA branches are clearly identified, while the TO and TA are 
barely visible due to a selection rule [3Bj. Note the very low 
intensity of the ZO mode close to T. Yellow lines are fits to 
the dispersion obtained with higher resolution data, where 
Kohn anomalies can be appreciated (see text for details). 



surements, respectively. Graphene was attained by dos¬ 
ing ethylene onto the clean Pt(lll) substrate held at 
1150 K. The completion of the first layer was reached 
upon an exposure of 3 • 1(F 8 mbar for ten minutes (24 
L, 1 L = 1.33 ■ 10 -6 mbars). The graphene layer was 
characterized by Raman spectroscopy, low-energy elec¬ 
tron diffraction, and scanning electron microscopy ex¬ 
periments I2Z1ESI- In particular, Raman measurements 
indicate the unique presence of monolayer graphene do¬ 
mains. 

HREELS experiments were performed by using an 
electron energy loss spectrometer (Delta 0.5, SPECS) 
with an energy resolution ranging from 1 to 3 meV. The 
dispersion of the peaks in the energy loss E\ oss was mea¬ 
sured by moving the analyzer while keeping the sample 
and the monochromator fixed. The phonon in-plane mo¬ 
mentum was determined from qu = ki sin 9i — k s sin 9 S 
as 

q\\ = y/2mE p ^sin0, - yjl - E\ oss /E p sin 0 S ^ (1) 

where 9i and 9 S are the incident and scattering angles, 
and we set h = 1. The impinging energy E p and the 


incident angle 0; were chosen so as to obtain the highest 
signal-to-noise ratio. A primary beam energy of E p =20 
eV provided the best compromise among surface sensi¬ 
tivity, the highest cross-section for mode excitation and 
momentum resolution. The angular acceptance of the 
apparatus was a = ±0.5°, which determines the momen¬ 
tum resolution as 

Aqy = sj2mE p ^cos0i + \Jl - E\ oss /E p cos 9 s ^j a (2) 

For the investigated range of qy, Aqy was found to range 
from 0.005 near T to 0.022 A -1 at K. To obtain the en¬ 
ergies of loss peaks, a polynomial background was sub¬ 
tracted from each spectrum. The resulting spectra were 
fitted by a Voigt line shape [28]. All measurements were 
made at room temperature. 

Experimental results - The phonon spectrum in 
graphene is composed by six phonon branches. Four vi¬ 
brate in-plane and are labeled transverse and longitudi¬ 
nal acoustic (TA and LA) and optical (TO and LO). The 
other two are acoustic and optical out-of-plane vibrations 
(ZA and ZO). Fig. [l] shows a map of the HREELS in¬ 
tensity as a function of energy and momentum, where 
four phonon branches (ZA,ZO,LA,LO) can be clearly 
observed. Their dispersion is consistent with previous 
works [25] [29] [30] ■ The other two phonon branches (TO 
and TA) are practically absent because a selection rule 
that forbids the emission of odd phonons under reflec¬ 
tions by the scattering plane 2B]. This is also clearly 
observed in Fig. [2j where the energy loss is plotted for 
a selected angle, and four peaks are clearly identified. 
The displacement pattern for each of these four modes is 



FIG. 2. Selected phonon spectrum from the dataset of Fig [T] 
at 9 S = 48°, which corresponds to q\\ = 0.56 A . The TA 
and TO modes are not observed due to a selection rule. Inset: 
Phonon displacements for the four phonons observed. Crosses 
and circles denote displacements in and out of the plane re¬ 
spectively. 
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FIG. 3. Phonon dispersions close to the V point, showing cusps around around q c ~ 0.13 A -1 . a) Dispersion of the LO phonon. 
A fit to Eq.|s](full line) gives an EPC A lo = 0.029, and parameters uj°lo = 200.6 meV, cllo = —110 meV A 2 , &lo = 558 meV 
A 4 , b) Dispersion of the ZO phonon. A fit to Eq. 0 (full line) gives an EPC Xzo = 0.087 and (jo° zo = 102.3 meV, azo = —810 
meV A 2 , bzo = 5243 meV A 4 , c) Representative spectra used to produce b). The frequencies are marked with a dashed line. 
Note the fits in a) and b) are also shown in Fig. [I] as yellow lines. 


shown as an inset to Fig. [2j 

We now focus on the dispersion of the ZO and LO 
modes around T. As usual, the ZO mode is significantly 
softened compared to the LO mode, due to the higher 
freedom for atom motion perpendicular to the plane with 
respect to the in-plane motion. The intensity of the ZO 
phonon at small momenta has very low weight, as shown 
in the intensity plot reported in Fig. 0 (see also [28]) . To 
enhance the signal-to-noise ratio from the ZO phonon in 
the nearness of T, each spectrum was acquired for several 
days. While this allowed us to resolve the ZO loss peak, 
the error bars for its frequency are still in general larger 
than for other branches. A careful fitting procedure |25] 
was used to extract the experimental frequencies for the 
LO and ZO phonons for small momenta, which are re¬ 
ported in Fig. |S3| The most striking feature in these 
plots is that they both display a clear cusp at the same 
momentum q ~ 0.13A . This strongly suggests that 

both cusps are Kohn anomalies at q = 2 kp, which origi¬ 
nate from the interaction of phonons with electrons. The 
Fermi wave-vector kp = Ep/vp can be estimated from 
ARPES measurements of graphene/Pt(lll) [23] . The 
reported values of the Fermi energy and Fermi veloc¬ 
ity are Ep ss 0.30 ± 0.15 eV and vp ~ 6 eVA. Thus, 
2 k F « 0.10T0.05A" , in good agreement with the posi¬ 
tion of the cusp found in this work. The position of the 
G peak in our Raman measurements [281 can be used to 
estimate ED 132] Ep ~ 0.42 eV, also in good agreement. 
While a Kohn anomaly for the LO mode is known to 
occur, the presence of the cusp for the ZO mode is unex¬ 
pected and represents the first evidence of the coupling 
of the ZO mode to electrons. In hindsight, by carefully 
inspecting the phonon dispersion recorded for graphene 
on Pt(lll) in a previous experiment by some of the au¬ 


thors ;29j , a dip for both LO and ZO at finite momen¬ 
tum is in fact seen but was not noticed or discussed. In 
Ref. [33] , the LO anomaly was incorrectly presented by 
shifting it to the F point due to a misunderstanding, but 
in fact occurs at finite momentum as well. To substanti¬ 
ate the claim that these cusps are in fact Kohn anomalies, 
we now compute the self-energies of both LO and ZO 
phonons following the conventions of Refs. Bang EH 
and compare the predictions to the experimental data. 

Theory - The dispersion of the phonons around T is 
modified because of their coupling to electrons, which at 
low energies can be modeled with a Dirac Hamiltonian 

H = J (vf°- k- /r) ip, (3) 

with <7 = (cr x ,cr y ) are Pauli matrices, vf is the Fermi 
velocity and /l is the chemical potential. This effective 
model is applicable up to energies Ap ss 1.5 eV or mo¬ 
menta A g = A e/vf = 0.25 A -1 , beyond which the dis¬ 
persion is no longer linear. Therefore, our predictions are 
only valid for phonon momenta within this range as well. 
The phonon Hamiltonian is 

H = J2j ■^2 U ite) b i,q b i, g ’ ( 4 ) 


with creation and destruction operators defined by the ef¬ 
fective displacements associated with each phonon mode 


Ui = 




where i = LO,TO, ZO and A c = 3v / 3a 2 /2 is the unit 
cell area, with a = 1.42A the nearest neighbor distance. 
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The dispersions of the phonons u>i(q) are analytic in the 
absence of EPC. For momenta close to the T point they 
can be expanded as w*(g) = w,° + aiq 2 + biq 4 , where ai®, 
cii and bi are parameters to be fitted from experimental 
data. The coupling between electrons and phonons is 
described by 

H e - P h = Fi j d 2 rui^Mii>, ( 6 ) 


with Fi the electron-phonon coupling. The matrix Mi for 
the different phonons is equal to Mlo = Q x <?, M to = 
q ■ a, Mzo = 03. Following Ref. [16], it is customary to 
introduce a dimensionless EPC as 

A ‘ - < 7 > 

The electron-phonon coupling induces a phonon self¬ 
energy Hi(q) that corrects the dispersion according to 

WRy = + cnq 2 + biq 4 + yll i(q/k F ). (8) 

In the static approximation, the self-energies for the dif¬ 
ferent phonon branches in terms of the dimensionless 
variable x = q/kp are given by [28] 


{^\Jl~ ^ + | arccos(2/:r)^ 9(2 - x), 

(9) 

nro(a;) = 0, (10) 

II zo(x) = ^ s ' v ^ (2 + £ arccos(2/a:)0(2 — x)) . (11) 


These expressions can be used to fit the cusps in the ex¬ 
perimental curves to obtain an estimate for the different 
electron-phonon couplings. To do so, we have fitted the 
dispersion parameters for both curves and determined 
the optimal Fermi level as Ep = vpq c /2 « 0.401 eV, 
which corresponds to a cusp mom entum of q c = 0.133 
A- 1 . The fits are shown in Fig. S3 with the experimental 
dispersion, and the result for the LO EPC is A lo = 0.029, 
in excellent agreement with the estimates obtained from 


Raman, \ LO = 0.027 — 0.034 [16]. The fit for the ZO 
gives A zo = 0.087, an even greater value. 

Discussion - Since this is the first work to observe the 
effects of electron-phonon coupling to the ZO mode in 
graphene, there are no measurements of A zo available 
to compare with. Theoretically, this coupling has only 
been estimated for a Si02 substrate in Ref. [19], giving 
a maximum value of F Z o = 7 eVA or A zo = 0.011 [55] , 
Our significantly higher value reflects the fact that the 
substrate is metallic and that, while direct hybridization 
with graphene is small, the surface electric field induces 
a stronger breaking of reflection symmetry which is re¬ 
sponsible for the EPC. 

Our finding that A zo has a significant magnitude im¬ 
plies that its effects should also be observable in other 


experiments. For example, the EPC is visible in the elec¬ 
tron dispersion, in the form of a kink at the phonon fre¬ 
quency that is observable by ARPES [9] . The existence of 
A zo would be responsible for an extra kink at lo Z o that 
could be resolved in future experiments. One could also 
expect a finite A zo to induce new Raman peaks. A peak 
at cozo from a first order Raman process (an analog of the 
G peak) for the ZO is not allowed because the ZO trans¬ 
forms as B 2l which is not contained in Ei x E\. However, 
a second order process (analogous to the 2D’ peak [IB]) 
is in principle allowed by symmetry. This process would 
give a peak at roughly 2oj Z Ch which unfortunately may 
be difficult to detect because of its overlap with the G 
peak. Therefore, either Raman or ARPES could provide 
an independent confirmation of electron-phonon coupling 
for the ZO mode. 

The finding that the substrate can induce an EPC to 
the ZO may also have important implications for trans¬ 
port in the high-field regime. The phenomenon of current 
saturation at high bias originates from inelastic scatter¬ 
ing with optical phonons, both intrinsic [36] or from the 
substrate [37] ■ Since the energy of the ZO (~ 100 meV) is 
significantly lower than the known phonons with strong 
EPC, the effects of the ZO may already be present in 
current experiments [3fiII3§| . and its presence should be 
accounted for in theory [40l - i42] . 

A related system where the standard (LO and TO) 
Kohn anomalies have been observed is graphene on 
Ir(lll) [3U]. According to ARPES [53], the Fermi level 
in this system is Ep ss 0.1 eV which corresponds to 
2k p sa 0.03 A . Kohn anomalies at such small mo¬ 
mentum are probably too difficult to resolve, which is 
consistent with the observed broadened dip around T in¬ 
stead of a cusp. Similarly, if the ZO EPC in this system is 
significant, such a feature should be observed for the ZO 
around T. The data in Ref. [30] are inconclusive on this 
matter. Finite momentum Kohn anomalies will be more 
likely to be found in graphene on metals more similar to 
Pt(lll), with weak hybridization but large charge trans¬ 
fer [441 . In metallic substrates with stronger hybridiza¬ 
tion such as as Ni(lll), Kohn anomalies are not present 
for the LO and TO phonons pTollTj] because the coupling 
of the graphene ir bands with the Ni d orbitals completely 
rearranges the electron bands. A Kohn anomaly for the 
ZO is therefore not expected in this type of system. 

In summary, in this work we have shown that the pres¬ 
ence of a substrate induces a significant electron-phonon 
coupling to ZO phonons in graphene, which is responsi¬ 
ble for a strong Kohn anomaly at q = 2kp. This finding 
paves the way to explore the many implications of the 
coupling of flexural phonons to Dirac electrons in sup¬ 
ported graphene samples and graphene-metal contacts. 
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I. REVIEW OF EXPERIMENTAL TECHNIQUES 


In this section we compare the different experimental techniques to measure phonons and explain why the ZO Kohn 
anomaly can only be observed with HREELS. There are four basic techniques to probe phonons: Neutron Inelastic 
Scattering, Raman scattering of X-Rays, Inelastic He Atom Scattering and Electron Energy Loss Spectroscopy (EELS). 

The first two techniques are better suited to the study of bulk 3D phonons. Neutrons are particularly well suited for 
bulk phonon dispersion studies because of their very small cross section with matter which allows them to penetrate 
deeply into the crystal. X-ray photons in grazing incidence, available at synchrotron radiation sources, are also a 
standard probe for bulk phonons and have in fact been used to study phonon dispersion in bulk graphite [7j. But 
none of this methods is suitable for studying phonon modes of a single layer of graphite because of the lack of surface 
sensitivity. 

Inelastic helium atom scattering is a more appropriate surface probe. However, in standard conditions the impinging 
energy is 65 meV. In order to observe modes at 100 meV it would be necessary to keep the nozzle at very high 
temperature, which implies a strong decrease of the signal. Inelastic helium atom scattering is more appropriate to 
study phonons up to 20-30 meV. The ZO phonon at around 100 meV has too high energy to be probed by inelastic 
helium atom scattering. 

HREELS measurements, on the other hand, have no limit in the energy loss and can easily access the high frequency 
range. They are surface sensitive measurements and have excellent resolution in both momentum and energy. This is 
thus the most suitable technique to study the full surface phonon spectrum, and in particular the only one capable of 
detecting the presence of the ZO Kohn anomaly. 
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FIG. SI. a) Selected HREELS spectra for different values of the parallel momentum transfer gy, evaluated by inserting the 
kinematic conditions and the loss energies in Eq. 1 of the main text. The intensity of the various phonons has been normalized 
to that of the elastic peak. Note the poor intensity of the ZO mode near T. b) EELS Intensity plot for phonon dispersion of 
graphene/Pt(lll), reproduced from Fig. 1 of the main text. Note again the very low intensity of the ZO mode close to T. The 
fits in Fig. 3 of the main text are displayed as yellow lines. 
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FIG. S2. Intensity of the ZO (filled squares) and the LO (empty circles) phonon modes as a function of the off-specular 
scattering angle, with the incidence angle fixed at 80° with respect to the sample normal. The intensity of phonon modes has 
been normalized to that of the elastic peak. 


II. FIT PROCEDURES 


In this section we present more details of the process employed to obtain the phonon dispersions displayed in the 
main text, and in particular in Fig. 3 of the main text where the Kohn anomalies are demonstrated. First, the 
intensity plot in Fig. 1 of the main text was produced from many different spectra acquired at fixed incidence angle 
(&i = 80°, all angles with respect to the sample normal) and at scattering angles varying from 9 S = 80° up to 15°. 
A representative subset of this series of spectra is shown in Fig. [Si] as a waterfall plot, along with the intensity plot 
which is reproduced here for ease of comparison. 

It is clear that finding the location of the Kohn anomalies in an intensity plot by simple inspection is not really 
possible due to the low intensity near T. To help the reader locate the Kohn anomalies in the dispersion of the LO 
and ZO in the intensity plot, we have also included the fits of Fig. 3 of the main text as yellow lines in Fig. |Sl[b). To 
quantify the low intensity, in Fig. [S2]we also present the integrated weight of the ZO and LO peaks as a function of 
the off-scattering angle with incidence angle kept at 80°. The intensity of the ZO mode increases with the off-specular 
angle, and thus with the momentum q \\, with a continuous behavior. Instead, the LO mode show maximum intensity 
for 15°-40° off-specular angles. Near T, i.e., near the specular conditions 9i = 6 S = 80°, the intensity of the LO and 
especially of the ZO mode is extremely poor. The spectra near T actually required several days of measurement in 
order to reach a sufficient signal-to-noise ratio. 

To quantify the dispersion of the LO and ZO around T accurately, we extracted the energy for each peak in each 


spectrum with a given gy with the following fitting procedure, exemplified in Fig. S3 First, an exponential background 


was subtracted from the raw HREELS spectrum. Then, the obtained spectrum was fitted with Voigt line-shapes. 
The Voigt line-shape is a convolution of Lorentzian and Gaussian line-shapes. The Lorentzian part takes into account 
the intrinsic line-shape of phonon modes, while the Gaussian line-shape is necessary to account for the experimental 
broadening due to the energy resolution of spectrometer. The extracted set of pairs (E, gy) was then fitted with the 
model dispersion relations discussed in the text. Both the data and the fits are presented in Fig. 3 of the main text. 


III. SAMPLE CHARACTERIZATION 
Scanning electron microscopy investigation 

In Fig. |S4| (a) we report a scanning electron microscopy (SEM) of the graphene/Pt(lll) sample, also reported 
elsewhere m■ The surface morphology appears to be homogeneous across all the sample (lxl cm 2 ). The SEM image 
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FIG. S3, a) Selected HREELS spectrum. The fit to an exponential background is reported with a dashed line, b) Fit procedures 
for TO+LO (left panel) and ZO modes (right panel) with the exponential background in a) substracted. 


shows a full coverage of the substrate surface by the graphene which forms a network of wrinkles (darker horizontal 
and vertical lines in Fig. [S4ja)). The wrinkles network develops on a micrometric length scale. Its morphology 
matches that obtained by low-energy electron microscopy measurements for graphene grown on Pt(lll) by carbon 
segregation from the Pt(lll) substrate and other metallic substrates [24] . 


Raman investigation 

The sample has been also analyzed ex-situ by Raman spectroscopy. Such experiments indicate that the surface 
of the Pt(lll) substrate is well coated with monolayer graphene, which is more ordered and homogeneous than for 
similar preparations previously reported in literature. A compressive strain is found for the graphene as the effect of 
the growth process. For more details on Raman investigation, see Ref. m- In Fig. |S4[b) we report a comparison 
between the Raman signal from the G peak in graphite and graphene/Pt(lll). Graphite shows a single peak (data 
taken from Ref. m ), which can be fitted by a single Voigt line-shape. By contrast, graphene/Pt(lll) shows a clearly 
asymmetric G peak, located at w = 1603 cm -1 . The position of the G peak can be used to estimate the Fermi 
level. A peak at 1603 cm -1 corresponds to a hole concentration [ 32 ] of n = 1.5 • 10 13 cm -2 , which corresponds to 
H = up’(?r7r) 1 / 2 ~ 0.42 eV, in very good agreement with the estimate obtained from the position of the cusp, and 
within the error of the ARPES estimate fj, = 0.3 ± 0.15 eV. 

It is interesting to note that the asymmetry of the Raman G peak could be originated in the emergence of a second 
order process with a peak around 2 u>zo, activated by the linear electron-phonon coupling. However, since this type 
of second order peak can probe a range of energies because momentum is not restricted to k = 0, further theoretical 
modeling is required to determine unambiguously that this peak indeed originates from the ZO. 


IV. THEORY: PHONON SELF-ENERGY 


All calculations are performed at T = 0. Room temperature (300 K) corresponds to 25 meV, which is much smaller 
than Ep or the optical phonon frequencies. With the definitions given in Eqs. EE in the main text, the phonon 
propagator is written as 


Gi(u>,q) = 


2 uji(q) 


to 2 - tOi(q) 2 - 2u} l (q)E l (uj,q) ’ 


(SI) 


where E,;(aj,g) = -ifHi{uj,q) is the self-energy, and i = LO,TO, ZO as in the text. The dispersion relation for the 


phonon can be obtained by solving for the pole in Eq. (SI) for small Aj, so that the renormalized phonon dispersion 
relation is given by 

UR,i{q) ~ tOi(q) + y H i(ui, q). 


(S2) 
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Wavenumber [cm ] 


FIG. S4. a) SEM image of a region of the graphene layer on a Pt(lll) substrate. Taken from Ref. [27]• b) Raman spectrum 
of graphite (data taken from Ref. [JTj) and MLG/Pt(lll) (data taken from Ref. 571 1 


which yields Eq. [8] in the text. The response function 11^ for a particular phonon is given by 

dcocPk 


Ki(E,q) = g SiV i 


(2tt) 


3 tr k)MiG(u) + E, k + q), 


(S3) 


where g StV = 4 accounts for spin and valley degeneracy, A4,; is the electron-phonon coupling matrix as defined in the 
text, and 


G(w,k) = 


co + g + a • k 


(w(l + ie) + g) 2 — k 2 ’ 

is the electron propagator, where we set vp = 1. Evaluating the trace in Eq. |S3| we obtain 

dcocPk Ni 


H-i(E, q) g SiV i 


(2tt) 3 [( w (l + ie) + /x)2 - k 2 ][((co + E)( 1 + ie) + g) 2 - |k + #] ’ 


where the numerators Ni for the different phonons are given by 

Nlo = 2 [ cj ( c o + E) — k • (k + q) + 2(/c • q)^ 1 / q 2 + 2k • q\, 

N T o =2 [co(co + E) — k ■ (k + q) + 2 (jc x q) 2 /q 2 ], 

Nzo = 2 [ui(u> + E) - k ■ (k + q)], 


The frequency integral is performed with the residue method, giving 


n i(E, q) — g StV 


d 2 k f 0(k-fi)Ni | g(|fc + 9l-^U- fi+ |t + ,1-, 

(2tt) 2 2\k\[(k + E(1 +ie)) 2 — + g| 2 ] 2|fc + q \[(|fc + g| + E(1 + ie)) 2 - \k\ 2 ] 


(S4) 


(S5) 


(56) 

(57) 

(58) 


(S9) 


Since these integrals are ultraviolet divergent as power counting reveals, it is customary to split them into the 
contribution coming from the undoped Dirac cone plus a finite correction induced by finite chemical potential. This 
is done by replacing 0(k — g) = 1 — 9(g — k) in the integrals, giving 

I ME, q) = n r°(E, q) - An i{E, q), 


(S10) 
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FIG. S5. Left: LOLo(q) f° r different values of g. Right: u>zo(o) for different values of /r. In both cases the coeficients of the 
analytic part of the dispersion are those obtained from the experimental fits in the main text. 


with 


All i(E,q) = g SjV 


d 2 k 


9(n - k)Ni 


0(/i \k + q])N i \ ui= _ E+ ^ + ^_^ 


2|fe|[(fc + E(l + ie)) 2 - \k + q\ 2 ] 2\k + q\[(\k + q\ + E( 1 + ie)) 2 - \k\ 2 ] 


nf ( E , q) can then be computed with dimensional regularization methods giving |5j, i l6] 

n To = 9 -^V~E 2 +q 2 , 


TTM-0 _ _ 
il TO — 


3s, i 


E 2 


16 ,J-E 2 + q 2 ' 


n^° =^fy/~E 2 + q 2 . 


(511) 

(512) 

(513) 

(514) 


The correction terms become complicated expressions but can be computed analytically in the static approximation 
(E = 0) giving [TU] 


An LO = 


47T 


qn 


— ( ^ arccos(^) + —\/q 2 — 4/x 2 ) 9{q — 2/z) 


aiito =o, 


An zo = 9s ’ 


qir 


4n L2 -2n-q&rccos( 2 f)9{q-2n) 


Combining the two contributions, Eq. |S10[ we finally get 


^ arccos( ^) + — \/q 2 — 4 [i 2 


iWE = o)=^ 

n to(e = o) =o, 

n zo (E = 0) 2/z + garccos(^)3(q - 2/x) 

47T L H 


9{q - 2/r), 


(515) 

(516) 

(517) 


(518) 

(519) 

(520) 


The results in the main text are obtained by recovering vf and setting x = q/irtFlj)). Note that the functional 
dependence of the /i = 0 part for the ZO was anticipated in Ref. |21j , and by symmetry arguments |l5j , it is the same 
as that of the TO phonon at the K point [IS] . 

These self-energies can be then fed into Eq.[S2]to determine the phonon energies. As an example, in Fig.[S5]we have 
plotted u>lo and wzo for different values of g,, with the coeficients of the analytic part as obtained from the fits to 
the experiment. This plots emphasize the dependence of the cusp position on the chemical potential, or equivalently, 
on kF- 


V. LIST OF ACRONYMS USED 


• EPC - Electron-phonon coupling 



























ZA - Out-of-plane acoustical 
TA - Transverse acoustical 
LA - Longitudinal acoustical 
ZO - Out-of-plane optical 
TO - Transverse optical 
LO - Longitudinal optical 

HREELS - High-resolution electron energy loss spectroscopy 
ARPES - Angle-resolved photoemission spectroscopy 



